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A number of versions of the extended three-band Hubbard models have been studied re- 
cently [|1|-|^ with the expectation that they are relevant to the physics of the high-T^ cuprates. 
The complexity of these models leads to a broad spectrum of phenomena, from metal-to- 
insulator transitions, charge and spin-density waves and superconductivity to more mundane 
single particle band-structure effects. Most of the studies searched for novel low-energy ex- 
citations which would explain the departures from Landau's Fermi Liquid phenomenology 
implicit in many of the experimental results in the normal state of the cuprates. 

In the simplest discussions based on mean field theories of strongly correlated Fermi 
systems such non-Fermi liquid behaviour occurs above some small energy scale, Tcoh- Already 
within the three-band models there are a number of different physical mechanisms resulting 
in the appearance of small energy scales, namely the proximity of the chemical potential to a 
van Hove singularity in the single-particle density of states , heavy Fermion behaviour |^ 
and the Brinkman-Rice metal-to-insulator transition and the effects of strong anti- 

ferromagnetic spin fluctuations in the paramagnetic state 

In this letter we focus on a band-crossing which arises upon hole doping (no > 1) for 
values of the — hopping matrix elements larger than a filling-dependent critical value, 
tpp > 0. In the vicinity of tpp = tpp the lowest (Cw-like) band is only weakly dispersing 
along the Brillouin zone boundary leading to a nearly one dimensional "extended" van Hove 
singularity in the density of states (cf. Ref. 0). In addition, once the doping is increased 
to a particular concentration, no = ficr, the band-crossing occurs at the Fermi level. (This 
rare situation in the context of three dimensional band-metals is best exemplified by the self- 
intersecting Fermi surface of graphite |10|.) Moreover, at n^, the resulting Fermi surface 
(FS) is perfectly nested with an incommensurate nesting wave vector. A special feature of 
this situation is that the nearly nested character of the Fermi surface survives for a wide 
range of dopings around ricr- In particular, in the vicinity of ricr the resulting low-energy 
scale varies with doping slower (Tcoh oc {n.Q — ncr)'^) than that expected from the saddle-point 
van Hove scenario (oc — |no — |/ In |no — ncr\)- 

We emphasize that in the experimentally interesting case of hole doping, the band- 
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crossing only occurs for a particular sign of tpp(> 0). Under these conditions the van Hove 
singularity of the density of states does not scan the Fermi level with increasing doping 
and thus the van Hove mechanism discussed in cannot be realized. Since, generically, 
the crossing of two eigenvalues of a hermitian operator requires the fine tuning of three 



real parameters |]TT| band-crossing in two-dimensions is typically avoided. Our situation is 
special as the mean field Hamiltonian can be reduced to a real matrix and thus crossing 
may occur through the tuning of only two parameters (e.g., the two components of the 
momentum, (k^jky)). We note that avoided crossing along a line in 2-dimensions (or at a 
point in 1-D) is responsible for the heavy Fermion regime of our mean field solution [|l^. 

In the present paper, we concentrate on the model with tpp > 0. We note that this 
choice disagrees with the situation which apparently takes place in all the cuprates, where 
the estimates extracted from LDA calculations [ll3| , p!^ imply that tpp is negative. In YBCO 



near optimal doping the corresponding FS agrees qualitatively with the ARPES results [jT5 
Some questions however arose in the context of BiSCO where two sheets appear in one of 
the earlier interpretations of the photo-emission results |T^ and there remains controversy 
of whether the origin of the two sheets is a lattice superstructure effect, bilayer splitting or 
some other effect. The main reason for studying a model with tpp > is the fact that, as we 
will see shortly, a rich array of many-body and single-particle phenomena arise in a natural 
way. It is interesting, however, that in some cases one is able to draw a parallel between our 
results and the experimental observations. If this model has any relevance it can only be 
so due to strong correlations which are believed to be crucial to the physics of the cuprates 



1% and are not properly accounted for within the LDA. To exemplify this point, imagine 
starting with a tpp = tight binding model with a large value of the Hubbard U and deriving 
an effective low-energy Hamiltonian by "integrating out" charge fluctuations on the Cu-sites 
one obtains an effective — hopping amplitude which changes sign with increasing U. In 
the ?7 ^ oo limit this yields t%f > (cf. Ref. 0). 

Our starting point is the extended three-band Anderson-Hubbard Hamiltonian describing 
the interaction between the oxygen Px,Py {Pxa,Pya) and copper dx2-y2 (da) orbitals. As 



depicted in Fig. 1, the Hamiltonian includes the direct hopping of electrons between oxygen 
sites (tpp) and a hybridization between copper and oxygen orbitals (tpd). Apart from the 
Hubbard repulsion at the copper sites - here after taken as infinite - we also take into 
account the nearest neighbor Coulomb interaction, V, between copper and oxygen. The 
infinite Hubbard repulsion is incorporated through the replacement of the original copper 
orbitals operators, dcr, by the projected Fermion operators, d^r = dail—n^cr) which eliminate 
double occupancy at the copper sites; here ria' = d\,dai is the local number operator for the 
orbital with spin a' . Below we treat this Hamiltonian by a simple mean field (Hartree-Fock) 
approximation: as explained in Refs. [^p!9|, the factorization of the equations of motion 
requires a generalization of the usual Wick's theorem. The resulting mean field equations 
can be adequately described by the quasi-particle Hamiltonian, 

Hmf = -2^t;, E dUk) [sin ^P.^k) - sin ^Py^k)] - 4tp, E PlMPyM sin y sin ^ + h.c. 
He', - /i) E di{k)d.{k) + - /i) Y.[PlAk)PUk) + Pl{k)Py.m . (1) 

k,a k,a 

Here t'^, = — n^itpd + V\/Stpd), = + A + AupV and e'p = ep + 2ndV are renormahzed 
values of the hybridization amplitude, and local energy levels of the copper and oxygen 
orbitals, tp,, e^, ep, respectively; given A*^°) = Cp — and the total filling, uq = + 2np, 
the values of the local copper (n^) and oxygen (rip) occupancy, the energy shift. A, and the 
chemical potential, n, are to be determined self-consistently from the mean field equations 
(MFE), 
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2^J^ / .A (e.(fc)-A)'-16t^pSin^^sin^» 
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Uj^i {(^iik) - (^j{k) 



, 8 / 2 ^ 3 ,^ 4tppSin^^ + Q(fc)-A . , ... 

where /i = /i — e'^, A = e^ — e^, N is the number of copper sites, and npix) = [exp 1]"^ is 
the Fermi-Dirac distribution function. Finally, the quasi-particle energies, uji{k) = e', + ei{k), 



3, should be determined from the cubic secular equation for ei{k), namely, 



2Ae^ + 



A - 4tp^(sm — + sm ^)-16tppSm — sm 



2 ~"' 2 

+4A t;^,(sin^ ^ + sin^ ^) - 32t'^,t,, sin^ |i sin^ | = 0. (5) 
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Equations ( ^ |D define the filling-dependent bandstructure. The no-double occupancy 
constraint manifests itself through two important effects: (i) the upward shift of the copper- 
like band and (ii) its associated band narrowing. The first occurs as a result of the no 
double occupancy constraint which forces the copper component of the lower band to lie 
within a hybridization width of the chemical potential; while the band narrowing is due to 



the v^l — Ud factor multiplying tpd and reflects the fact that coherent hybridization involves 
rare charge fluctuations at the copper sites. The interplay between these correlation effects 
and the band structure leads to the rich phase diagram shown in Figs. 2. To summarize 
the possible behaviours imagine starting with the bare copper level below the oxygen level, 
A(°) > and increasing the filling from no = 0. We will consider the A*^°^ vs. uq plane for 
two qualitatively different cases, tpp < to and tpp > to, where to = \^cn\tpd\ + '^2^ with the 
two constants given by ai = (4-\/2 — 5)/it ^ .21 and 02 = [vr — 4(v^ — l)]/47r ^ .12. In 
each case we single out three A'^°^ = const cuts, denoted by A, B, and C in the figures. For 
tpp < to and for small values of A'^^-' (cut A in Fig. 2a) one starts with conventional metallic 
behaviour for uq << 1. Below half filling one then scans the the chemical potential through 
the van Hove singularity (saddle point) switching from hole-like to electron-like FS. The 
situation at no < 1 is reminiscent of what would happen at no > 1 if we reverse the sign of 
the oxygen-oxygen hopping. The low energy scale is provided by the difference between the 
Fermi level and the saddle point energy (see fig. 3), and although in the case of the Fermi 
energy equal to e^Hs one does get the linear temperature dependence of the quasiparticle 
relaxation rate [0, the corresponding crossover temperature T^, above which the relaxation 



rate becomes linear in T is in this case expected to scale as (no — riyHs)/ In |^o — n^Hs\ pO 
where n^us is the value of filling corresponding to fi = e^us- The physical consequences of 
the presence of a small energy scale fl — e^us have been discussed in the literature in great 



detail M\. 



With further increasing no the chemical potential scans through another van Hove sin- 
gularity associated with the bottom of the lower oxygen-like band, which, at T = 0, causes a 
negative jump in the compressibility, dno/d^. Above this doping concentration, drio/dfi < 0, 
and the system is thermodynamically unstable with respect to phase separation. Within the 
region of instability, the absolute value of the (negative) compressibility increases and, at 
some point, diverges. The compressibility then changes sign and the system enters a normal 
metallic phase with large and rapidly decreasing value of the compressibility. 

The shape of the unstable region is highly sensitive with respect to the parameters of 
the Hamiltonian. The phase separation would be circumvented altogether in the presence of 
a real long-ranged Coulomb interaction. Beyond half-filling the bottoms of the two lowest 
bands coincide ("tangency") beyond which a crossing of the Cu- and lower 0-like bands 
(conic point) emerges. At sufficiently large filling the chemical potential scans through the 
band crossing energy. Along cut B in Fig. 2a a new feature arises at tiq = 1 (between the 
two van Hove singularities), namely the Brinkman-Rice metal-to-insulator transition (BRT) 
where the copper-like band becomes completely fiat (see below). Finally, along cut C, in 
addition to the BRT, a heavy fermion regime with exponentially small 1 — and a conic 
band crossing point below the chemical potential appears immediately above half filling, 
extends over a finite range of dopings and then terminates via a smooth crossover into a 
conventional metal. In the case of cut C, the phase separation region is located above this 
crossover. In the metallic phase, the conic point eventually crosses the chemical potential. 
For tpp > to both van Hove singularities and the tangency appear already below half filling. 
The only other qualitative difference from the tpp < to case is the possible appearance of 
reentrant heavy fermion behaviour below the BRT critical value of A'^^^ (see cut B in Fig. 
2b). Numerical calculations show that the features of the phase diagram persist at finite 
temperature, with the value of A'-''-' decreasing as temperature increases. We are now in 
position to discuss these features in more detail. 

The mean-field bandstructure is characterized by three energy scales, namely the band- 
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width of the lower band, the difference jl — e^Hs between the Fermi level and the van Hove 
singularity (corresponding to the saddle point in the lowest band dispersion) and the dif- 
ference between the Fermi level and the energy ecr = —t'pd/tpp of the band crossing point. 
The values of these three quantities, computed along the three cuts A, B, and C of fig. 2a 
are plotted in fig. 3. One can see, that indeed the lowest energy scale in the metallic phase 
above half filling over an extended range of dopings is provided by /i — tcr, which vanishes 
at some particular filling ncr(A^°-') and at small values of |no — ncr\ ^ 1 depends linearly on 
filling, /i — ecr oc — no. 



The FSs that emerge when the Fermi level lies close to are shown in Fig. 4 . It is a 
striking feature of the present model that when the band crossing occurs exactly at the Fermi 
level, the FS is formed by the straight lines = k^r and ky = k^r (with k^r = — T^ncr/^)-, 



and is therefore perfectly nested along the coordinate axes. It is well-known [|2^ that in 
such a situation the quasiparticle relaxation rate becomes linear in temperature. One can 
therefore expect that as the value of filling approaches ricr-, the quasiparticle relaxation rate 
undergoes a crossover to the linear temperature dependence; this is also the case when, at 
some fixed value of no, temperature increases beyond certain crossover value T^. 

The simplest way to exemplify this behaviour is to replace the nested parts of the FS by 
two nearly straight segments of length L and curvature radius L, also assuming that the 
velocities of quasiparticles on these segments are equal to v and antiparallel to each other. 
One can then estimate the contribution of these segments to the off-shell decay rate due to 
a weak interparticle contact interaction, U5{r — r '). For uo > vL'^/R, and when momentum 
vector p lies at the middle of the straight segment of the FS, the decay rate is linear in 
frequency, ImS(u;,p) ~ (SU'^L'^uj)/ {Wn^v"^) whereas at very low frequencies, u <^ vL'^/R, 
one gets 

RU' 



ImS(ti;,p) 



ojR 3 



(6) 
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The finite-temperature on-shell quasiparticle relaxation rate can be estimated by substitut- 
ing T ^ uo m these expressions. The quantity L^/ R measures the deviation of the nearly-flat 



segments of the Fermi surface from the tangent straight hnes and provides a new low-energy 
scale, vl?lR. As expected for values of oj.T » vL'^/R ImS is proportional to 

X = max{u;(> 0), T}, whereas at low energies the behaviour reduces to the 2D Fermi Liquid 
result ImS oc — a;^lna; ||24|| . 



For our mean field band structure R oc {fl — ecr)'"^ in the limit of \jl — tcrlli\ << 1 and 
L ~ TT. Since jl — tcr oc ricr — no, it follows that the crossover scale T* = vL'^/R oc (no — ncr)^. 
As a result, the linearity of the relaxation rate survives over a wider range of dopings than 
in the van Hove scenario for which oc {{uq — ricr) /ln{nQ — nc,)\. 

Let us now take a closer look of the band crossing point. Expanding the dispersion law 
about this point leads to the following two branches. 



ei,2(fc) = A{k) T \jB{k){K + K- 2K,y - C{k){K ~ kyf , (7) 

— * 

where A, B and C are smooth functions of momentum k. Note that the square-root be- 
haviour in (|^) is associated with the non-analytic behaviour of the matrix elements of the 
unitary transformation which diagonalizes the mean field Hamiltonian; more precisely, the 
limiting values of these coefficients aX k ^ k^r depend on the direction of approach. This af- 
fects, for example, the interband matrix element of the quasi-particle position operator which 
diverges at the crossing point as ri2{k) oc (ei(A;) — e2{k))~^. On the other hand, the corre- 
sponding interband matrix element of the velocity operator Vi2 = {df/dt)i2 = i{ei — e2)ri2 
remains finite at A; ^ kcr but its limiting values still depend on the direction of approach. 
In turn, the leading interband contribution to the optical conductivity (ignoring all quasi- 
particle interactions) [p5[] . 



ReaS = ^ [ vUkHiik) \np{eS)) - nF{e2{k))\ 5{e^{k) - e2{k) + uj)(fk , (8) 



□ 2 

(the indices a, (3 take values x or y) displays an anomaly associated with the peculiar be- 
haviour of Vi2 in the presence of band crossing. The integration in ( ^ leads to the result 
sketched in fig. 5. Note the presence of two square-root singularities in the frequency 
derivative of the conductivity: one at the threshold frequency 

8 



r, ^pd + tpp^ I I 

below which RecT^|(co') is equal to zero, and another at = Ut ■ (St^^ + 2tppA)/tp^. These 
singularities originate from the tangency between the ellipse 62 (^) = ei(^) + ^ and the two 
borders of the stripe fi < ei{k) < fl + u . At larger frequencies, u ^ ut, uj[ (but u still much 
smaller than the bandwidths), 'Rea^'^iu) approaches a constant value. Finally, if the band 
band crossing occurs precisely at the Fermi level, the independence of Re(Ji2(co') on frequency 
survives down to = 0. It is amusing to note the resemblance of the threshold behaviour 
in Fig. 5 with the "mid-infrared peak" observed in optical conductivity experiments ||26|| . 
Within this scenario the weak doping dependence of the "peak" position could only be 
accounted for if the system is sufficiently far below the critical filling, Ucr- 

It is easy to see that the value of the threshold frequency and the nature of singularities 
at a; = uJt.,uj[ is in fact unaffected by the RPA corrections. The latter, in principle, might 
add excitonic features to the profile of Heaapiyj) (delta-functional peaks below the threshold 
or Lorentzian peaks above the threshold). We have checked, however, that excitonic bound 
states do not occur for physically relevant interaction strengths (V^ < W'^ /ut-, where W is 
the conduction electron bandwidth). Also, the effects of resonances above threshold is small 
for uj ojt- Thus, at frequencies of the order of ujt RPA modifies the rigid-band result by a 
factor of order unity. 

We close by noting that when the Fermi level lies in the vicinity of a band crossing 
point, one should be able to observe the phenomenon of magnetic breakdown [^], namely, 
the restructuring of the FS that takes place at sufficiently high magnetic fields through the 
tunneling of electrons between different FS sheets. This phenomenon and its consequences 
for magneto-transport in weak fields in the presence of a conic point will be described 



elsewhere 28 
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FIGURES 

FIG. 1. The arrangement of copper d-orbitals and oxygen p-orbitals in the Cu02 plane, with 

signs "+" and accounting for the pliasc factors in tlic atomic wave functions, wliich result in 
the alternating signs of hopping terms in the tight-binding Hamiltonian. 

FIG. 2. The phase diagrams for tpd = 1.3 eV, tpp = 0.65 eV, and V = 1.25 eV (a), and for 
tpd = 1-3 eV, tpp = 0.8 eV, and V = 0.5 eV (b). Notice the difference in the shapes of the 
heavy fermion (HF) region. The solid bold lines represent BRT, the dashed bold lines - a smooth 
crossover between HF and normal metal behaviour. The dotted line in fig. 3b corresponds to 
analytical result for the HF to normal metal crossover, obtained in the limit of small no — 1. 

FIG. 3. The bandstructure energy scales: the bandwidth of the lowest band (solid lines), jl — Ccr- 
(dashed lines), and jl — CyHs (dotted lines), computed along the three cuts A, B, and C in fig. 2a. 

FIG. 4. The Fermi surfaces for no = 1.23 (dotted lines), no = 1.46 (no ~ Ucr, solid lines), and 
no = 1.64 (dashed lines). The bare parameters of the system are tpd = 1.3 eV, tpp = 0.65 eV, 
y = 1.25 eV , and A^^) = 3.26 eV (cut B in fig. 2a), and the number 0, 1 and 2 indicate the 
number of filled bands in each region of the Brillouin zone. 

FIG. 5. The rigid-band estimate for the interband term in the optical (a6)-plane conductivity 
(ohm~^cm~^) versus frequency (eV) at no = 1.44, tpd = 1.3 eV, tpp = 0.65 eV, V = 1.25 eV , and 
A(°) = 3.26 eV. We used the value of 12 A for the interplane spacing. 
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